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1. Introduction 

Study of fluctuation and correlation have special importance while investigating the existence 
of possible phase transition in relativistic heavy ion collisions. Event by event fluctuations in the 
thermodynamic quantities provide important insight towards the physical properties of the dense 
nuclear matter produced in the collisions along with the knowledge of nature of quark-hadron phase 
transition The enhanced fluctuation of energy density points to first order phase transition 
whereas a second order phase transition may result into divergence of specific heat. Study of event 
by event fluctuation in (pr), charged and neutral multiplicity correlation are tools those could be 
used to look for the nature and order of the phase transition. Another aspect of these studies is to 
search for the existence and location of the critical point on the quark-hadron phase diagram [Q]. It 
is believed to observe large fluctuations of thermodynamic quantities in the vicinity of the critical 
point. 

Now, the measured fluctuation from the relativistic heavy ion collisions may have contribu- 
tions from many different sources. Firstly there are quantum fluctuations, which arise if the specific 
observable does not commute with the Hamiltonian of the system under consideration. But these 
fluctuations probably play a less significant role in heavy ion collisions. Secondly and most im- 
portant in our case are the so called "dynamical fluctuations" reflecting the dynamics and response 
of the system under external forces. They help to characterize the bulk properties of the system 
like the instability arising at the phase boundary if the system undergoes a phase transition. Finally 
there are contributions from known sources including (a) geometrical (impact parameter, partici- 
pant number, detector acceptance), (b) energy momentum and charge conservation, (c) anisotropic 
flow, (d) Bose-Einstein correlations, (e) resonance decays, (f) correlations from jets and minijets 
etc. These contributions are needed to be understood and subtracted in order to estimate the true 
extent of dynamical fluctuation that tell us about the properties of the system. One or more of the 
following methods are used to achieve this 

• construction of mixed events which contain the known fluctuation contributions 

• simulation of known fluctuations using the physical phenomena those contribute to them 

• formulating proper quantity which is robust against contributions from known sources 

With the advent of the recent experiments with relativistic heavy ion collisions many fascinat- 
ing results are available for the understanding of the QGP phase transition in particular as well as 
the overall nature of the dense matter in general. In this review the results obtained from Solenoidal 
Tracker At RHIC (STAR) experiment will be discussed. 

2. STAR experiment 

STAR experiment comprises of a Time Projection Chamber (TPC), as its central tracking 
detector, with a pseudorapidity coverage of \r\ \ <1.0 and azimuthal coverage —n<<p<n. All 
results discussed here in this review are obtained from analysis of data taken by the TPC. The 
charged tracks are identified by using their relative energy loss in TPC. Collision centralities are 
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measured by the number of charged particles with the pseudorapidity range |?]| < 0.5. Glauber 
Monte Carlo has been used to estimate the number of participants involved in the collisions for 
each class of centrality. 

The data used for the results in this review are from Au+Au collisions at y/s^ = 20, 62.4, 130 
and 200 GeV. 

3. (pt) correlation and fluctuation 

(pt) of the particle produced in a collision corresponds to the temperature of the system. 
So the event by event fluctuation in (pj) were proposed to estimate the fluctuation in tempera- 
ture. But recent (px) fluctuation studies reveal contributions from many other sources including 
flow and minijiets to it. In STAR there are several analyses studying the event by event fluctua- 
tion/correlation of (pt). 



Figure 1: (pj) fluctuation. Upper panel: 
Event frequency distribution of (pj) fluctu- 
ation measure from data (histogram) com- 
pared to statistical references, such as, F dis- 
tribution (dashed curve), Monte Carlo refer- 
ence (solid curve underlying F distribution). 
The broadened distribution (solid curve un- 
derlying data) is obtained by raising it to a 
power calculated from the measured value of 
fluctuation from statistical analysis (see text). 
Lower panel: Difference in upper panel, be- 
tween data and F reference (histogram) and 
difference between F reference and broad- 
ened distribution (solid curve). 

Fig. |l| shows a 14% deviation of event by event (pj) distribution from the statistical reference 
constituted (a) from F function distribution and (b) from Monte Carlo simulation ^ for central 
Au+Au collisions at =130 GeV. The broadened F function distribution is obtained by raising 
it to a power calculated from the measured value of fluctuation, Aa p , : „, from the statistical analysis. 
The agreement of this broadened distribution with the data histogram ensures the fact that the 
measured fluctuation excess is a property of most particles in most collisions, is not produced by a 
minority of outliers. The excess in width of data over the reference occurs rather smoothly on both 
sides of the peak. Absence of significant deviation in this rules out evidence of anomalous event 
classes, as might be expected from critical fluctuations. In Fig. || charge dependent and charge 
independent (pi) difference factors are plotted in dependence of event centrality. This shows that 
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the charge independent dynamical fluctuation is more than the charge independent ones and it 
varies smoothly and monotonically with collision centrality. 



Figure 2: Centrality dependence of (pj) fluc- 
tuation. Charge independent (solid triangular 
points) and charge independent (open triangu- 
lar points) difference. Parameterization (dashed 
curves) extrapolation of parameterizations to 
true primary particle number (solid curves), sys- 
tematic uncertainties (color bands). 

Another analysis in STAR [Q] uses two particle correlation approach as a measure of dynamical 
fluctuation. The (pj) correlation show non-zero values for all collision energies available at RHIC. 
Fig. |3| shows that for all collision energies the scaled value of correlation increases rapidly with 
centrality and then saturates for central collisions. To study the possible changes in the dynamical 
fluctuation with the change in the event average of (pj) (({pj}}) with centrality and incident energy, 
the dynamical fluctuation scaled by ((pt)) is plotted in Fig. |as a function of centrality as well 
as collision energies. Results from HIJING model calculation underpredict the measured values of 
correlation and do not predict the observed centrality dependence. 
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Figure 3: (dN /dt])(Ap t .iAp t .j) as a function of cen- 
trality and incident energy for Au+Au collisions com- 
pared with HIJING results. 
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Figure 4: ^/ {Ap,jAp t j) / ({pt)) as a function of 
centrality and incident energy for Au+Au colli- 
sions compared with HIJING results for corre- 
sponding systems. The inset shows the excita- 
tion function for the most central bin. 



There is another approach in STAR to study the scale dependence of (pj) correlations and 
associated angular autocorrelations [Qj. In this approach the scale dependence of (pj) fluctuations 
in (tj,0) is measured and the corresponding autocorrelations in difference variables (8ri,8<p) are 
obtained by inverting the (pj) distributions. These autocorrelations reveal interesting structures 
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in phase space, which are attributed to production of minijets in relativistic heavy ion collisions. 
Fig. IH shows the results from (pj) scale dependence for Au+Au collisions at y/SNN = 200 GeV. 



Figure 5: Distributions of 
Ap/7p^7(T7A,0 A ) for 80-90% 
(upper left), 45-55% (upper right), 
20-30% (lower left) and 0-5% 
(lower right) of total cross section. 
Monopole (constant), dipole and 
quadrupole components have been 
subtracted. 
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There has been a continuation of this approach by using the variance difference in (pj) angular 
correlation to determine the centrality and collision energy dependence of this quantity in STAR 
[f]]. Fig. H shows the results from this analysis. The left panel shows the centrality dependence 
of (pt) fluctuation for all collision energies at RHIC. The centrality has been measured as v ~ 
ZNbin/Npart where Nbin is average number of binary collisions and N part is the average number of 
participant nucleons. The vertical line at right estimates v for most central collisions (b=0). The 
lines are drawn to guide the eye. In the right panel the energy dependence of (pj) fluctuation for 
most central collisions is plotted. The full and open symbols correspond to data with and without 
SSC (HBT, Coulomb) correction. The solid curve is proportional to ln(y / S^v/10). 




Figure 6: Scale dependence of (pj) angular correlation, centrality dependence (left), energy dependence 
(right). 
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4. Net Charge fluctuation 

Fluctuation in conserved quantity like electric charge has been proposed to be a probe to the 
know about the QGP phase. There have been a number of theoretical investigations predicting the 
value of charge fluctuations for confined hadronic phase as well as deconfined QGP phase M. 
With the wide acceptance TPC of STAR the dynamic net charge fluctuations have been measured 
for Au+Au collisions for all available collisions energies. 

Fig. [7] shows the centrality dependence of the dynamical net charge fluctuation measured 
through a correlation function approach [||] for Au+Au collisions at y/stm =130 GeV JT0|]. There 
is a strong centrality dependence of the dynamical fluctuation and the two particle correlation might 
have been modified in central collisions relative to peripheral ones. The measured value of fluc- 
tuation is closer to the prediction for resonance gas. Fig. |8] shows that the dynamical net charge 



fluctuation is invariant of collision energies starting from top SPS energies [11]. 
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Figure 7: Centrality dependence of 
dynamical net charge fluctuation for 
Au+Au collisions at ^/sJw = 130 
GeV. 



Recent lattice calculations [|12[ |13|] predict an enhanced charge fluctuations in the hadronic 
phase where as suppression of the same at the high temperature phase of QGP. Moreover, these 
calculations observe prominent structures in the higher order moments of net charge distributions 
for temperature close to the transition temperature. Fig. ^ shows net charge distributions from 
charged particles with pi below 1 GeV/c for Au+Au collisions at = 200 GeV [ 14] for differ- 



ent centralities. Analysis of higher moments of net charge distribution is underway to extract new 
information. 



5. Particle ratio fluctuation 



Large difference in the relative production of kaons to pions has long been predicted due to the 



difference in free enthalpy of the hadronic and QGP phase [15]. Two of the most interesting and 
widely discussed experimental observations in the area of event by event fluctuation in relativistic 
heavy ion collisions are the strangeness fluctuation measured in terms of ratio of Kaons over Pions 



at SPS energies by NA49 Collaboration Qlofl as well as the collision energy dependence of the 



measured fluctuation [17] 



STAR experiment has also measured the event by event fluctuation in K/71 ratio in Au+Au 



collisions in all available collision energies at RHIC [18]. The kaons and pions are identified 
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Figure 8: Excitation function for dynamical net 
charge fluctuation. Data from SPS also plotted 
for better comparison. 
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Figure 9: Net charge distributions from charged 
particles with pj below 1 GeV/c for Au+Au col- 
lisions at v?snn = 200 GeV 



from the relative energy deposition by the charged tracks in TPC. The number of Kaons as well 
as pions are calculated event by event and the ratio is obtained. The dynamical fluctuation is 
extracted by subtracting the statistical reference, constructed by mixed event technique, from the 
measured value of fluctuation. Fig. |K] shows the excitation function of dynamical fluctuation. 
The data points for NA49 have also been shown in this plot for comparison. Fluctuation in K/k 
remains independent of the collision energies available at RHIC. The centrality dependence of 
dynamical fluctuation in Au+Au collisions at -y/SAw = 62.4 and 200 GeV are also studied using 
a correlation function approach similar to the net charge fluctuation analysis. Fig. 11 shows that 
the measured value of dynamical fluctuation decreases with increasing collision centrality. This 
dilution of fluctuation could be attributed towards increasing production of correlated pairs from 
resonances as the collision becomes more and more central. 



6. Balance function 

One interesting probe to know about the history of hadronization after the relativistic heavy 
ion collisions is balance function Jl9|]. In the collision when the particles and anti particles are 
produced in pairs, they are correlated in coordinate space. If the hadronization occurs early, the 
charge - anti charge pairs are expected to be separated in rapidity due to expansion of the fireball as 
well as the re-scattering in the strongly interacting matter. On the other hand if the pair is produced 
in later time the the correlation is stronger. As the width of the correlation reflects the time of 
hadronization, this could signal to creation of high temperature QGP phase. 

The width of the balance function for unidentified charged particles as well as identified pions 



(Fig. 12) in Au+Au collisions at Wsnn =130 GeV are observed to be narrower in case of more 



central collisions [20]. In spite of sensitivity of this observation to effects like flow, resonance pro- 
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Figure 10: Excitation function for the dynamical fluc- 
tuation in Kj it ratio. The data from SPS are plotted 
for comparison. 
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Figure 11: Centrality dependence of dynamical 
fluctuation. 



duction and diffusion, it also indicates towards possibility of late hadronization when the collision 
is more violent. The width predicted from HIJING model calculation is consistent only with results 
obtained in peripheral collisions. 
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Figure 12: Width of the balance function for unidentified charged particles (left) and identified pions (right). 
The result from HIJING calculation in both cases are shown in the shaded band. 



7. Forward backward multiplicity correlation 

Short and long range (in rapidity) multiplicity correlations are predicted as a signature of string 
fusion. From the relationship between the multiplicities in the forward and backward pseudorapid- 
ity regions, a correlation coefficient could be defined as, 

b _ <N f N h >-<N f ><N b > _D 2 bf 
<Nj>-<N f > 2 D 2 ff 
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where Nf and Nb are the multiplicity in the forward and backward region respectively and D^, 
Dyy are the backward-forward and forward-forward dispersions. When strings fuse a reduction of 
the long range forward backward correlation is expected. Fig. |l3| shows the centrality dependence 
of backward-forward and forward-forward dispersions measured in Au+Au collisions at y^\w = 
200 GeV. These results are obtained by considering the charged particles within 0.8 < |rj| < 1.0 
[22]. The measured values are also compared to those from Parton String Model (PSM) calculation 
with and without string fusion. The reduction in the long range multiplicity correlation, D^, even 
compared to the collective PSM, is a possible indication that additional collective effect is present 
in central Au+Au collision at Ws^m = 200 GeV. 
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Figure 13: D^t, D 2 ^ as a function of centrality for Au+Au collisions at ^s^n = 200 GeV. Data are compared 
with PSM calculations with and without string fusion. 



8. Summary and outlook 

The wide acceptance of the STAR detector both in pseudorapidity (|rj| < 1.0) and azimuth 
(—71 < <p < 7i) allowed us to study event by event fluctuations in various observables like (pr), 
net charge, K/n ratio as well as correlation analyses with balance function and forward-backward 
multiplicity correlation. Results from (pj) fluctuation analysis show substantial excess in fluctua- 
tion, over the statistical reference, which scales smoothly and monotonically with centrality. (pi) 
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analysis in finer (tj,0) sub-space shows interesting structures which could be attributed to produc- 
tion of minijets at RHIC. The centrality dependence of multiplicity scaled value of (pj) correlation 
may show signs of effects such as thermalization, the onset of jet suppression. Net charge fluctua- 
tion is greater than the predicted value from charge conservation limit but closer to resonance gas 
picture. Amount of dynamical net charge fluctuation remain essentially unaltered as a function of 
collision energies available at RHIC. Results from analysis of higher moments of net charge are 
promising but needs more work to obtain a more comprehensive understanding. STAR obtains the 
first results from event by event fluctuation in K/n ratio at RHIC. This fluctuation is positive for 
all collision energies available at RHIC, however it remains practically independent of the collision 
energy. The dynamical fluctuation in K/n ratio decreases as the collision become more and more 
central. Results from balance function analysis indicates towards a delayed hadronization after 
the heavy ion collisions at RHIC. A reduction in long range multiplicity correlation, compared 
to model calculation, suggests presence of additional collective phenomena in central heavy ion 
collisions. 

One of the most important aspect of the fluctuation studies in relativistic heavy ion collisions 
is to search for the location of the critical point on the quark-hadron phase diagram. This could 
be achieved by scanning the QCD phase diagram in terms of temperature and baryonic chemical 
potential. A program has recently been proposed to run RHIC at low energies to probe the region 
of low temperature and rather high baryonic chemical potential for this purpose. The STAR experi- 
ment has a shinning future in this program. With the addition of a full barrel Time Of Flight (TOF) 
detector, STAR will increase the momentum range of identification of K and % to play a key role in 
studying strangeness fluctuation near the phase boundary. It is also expected to have more insights 
from the energy dependence of balance function at lower energies. 

With the continued developments of new analysis techniques as well as upgrades of the de- 
tector subsystems, STAR will continue to provide more and more comprehensive knowledge about 
the quark hadron phase transition in relativistic heavy ion collision through correlation fluctuation 
studies. 

References 

[1] H. Heiselberg, Physics Reports 351 (2001) 161. 

[2] M. A. Stephanov, K. Rajagopal and E. Shuryak, Phys. Rev. Lett. 81 (1998) 4816. 
[3] J. Adams et al. for STAR Collaboration, Phys. Rev. C 71 (2005) 064906. 
[4] J. Adams et al. for STAR Collaboration, Phys. Rev. C 72 (2005) 044902. 
[5] J. Adams et al. for STAR Collaboration, J. Phys. G 32 (2006) L37. 
[6] J. Adams et al. for STAR Collaboration, nucl-ex/0605022 (2006). 
[7] M. Asakawa, U. Heinz and B. Miffler, Phys. Rev. Lett. 85 (2000) 2072. 
[8] S. Jeon and V. Koch, Phys. Rev. Lett. 85 (2000) 2076. 
[9] C. Pruneau S. Gavin and S. Voloshin Phys. Rev.C 66 (2002) 044904. 
[10] J. Adams et al. for STAR Collaboration, Phys. Rev. C 68 (2003) 044905. 



10 



Fluctuation studies in STAR 



Supriya Das (for the STAR Collaboration) 



[11] C. Pruneau for STAR Collaboration, Nucl. Phys. A 774 (2006) 651. 
[12] S. Ejiri, F. Karsch and K. Redlich, Phys. Lett. B 633 (2006) 275. 
[13] R.V. Gavai and S. Gupta, Phys. Rev. D 72 (2005) 054006. 

[14] T. Nayak, Proceedings of International Conference on Strangeness in Quark Matter (SQM06), 
nucl-ex/0608021 (2006). 

[15] J. Rafelski and B. Miffler, Phys. Rev. Lett. 48 (1982) 1066. 

[16] S.V. Afanasiev et al. for NA49 Collaboration Phys. Rev. Lett. 86 (2001) 1965. 

[17] G. Roland for NA49 Collaboration, J. Phys. G 30 (2004) S1381. 

[18] S. Das for STAR Collaboration, Proceedings of International Conference on Strangeness in Quark 
Matter (SQM06), J. Phys. G 32 (2006) S541. 

[19] S. A. Bass, P. Danielewicz and S. Pratt, Phys. Rev. Lett. 85 (2000) 2689. 

[20] J. Adams et al. for STAR Collaboration, Phys. Rev. Lett. 90 (2003) 172301. 

[21] GD. Westfall for STAR Collaboration, J. Phys. G 30 (2004) S345. 

[22] T.J. Tarnowsky for STAR Collaboration, Proceedings of 22 nd Winter Workshop in Nuclear Dynamics, 
nucl-ex/0606018(2006). 



11 



